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Abstract 

We compute the spontaneous staggered polarization of the cyclic SOS model at the 
^. \ thermodynamic limit. We use the determinant representation for finite-size form factors 

obtained from algebraic Bethe Ansatz. 
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1 Introduction 

In a previous paper [9], we have initiated the study of form factors and correlation functions 
of the cyclic solid-on-solid (CSOS) model [HQ2] by means of algebraic Bethe Ansatz (ABA). 
In particular, we have obtained determinant representations for the scalar products of Bethe 
states and the form factors in finite volume. The aim of this paper is to show how to apply 
these results to the computation of physical quantities at the thermodynamic limit. As a 
simple example, we consider here the spontaneous staggered polarizations. 

The CSOS model is a two-dimensional square lattice with interactions around faces (IRF). 
On each vertex on the lattice, there is a height s which may take L different values, so that 
heights on adjacent sites differ by ±1. The difference of heights between two adjacent sites 
of the lattice is hence described by a variable e = ±1 attached to the corresponding bond. 
In the cyclic SOS model [HI 112] . heights are periodic in L (i.e. there exists so E C such that 
s — sq E li/Lli). There are six different allowed configurations around a face, of the type 
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such that e; + e,- = e' + e'-. 



J = s + e • + e^- 

To these configurations correspond six statistical weights W( , s fl". 6 ?.,.) which can be pa- 
rameterized in terms of elliptic theta functions. They depend on a parameter 77 (crossing 
parameter) which, in the cyclic case, is a rational number: 77 = r/L with r, L coprime inte- 
gers. Hence, the statistical weights of the CSOS model are L-periodic functions of the height 
s. 

It has been shown in |12tlll| that the transfer matrix of the model with periodic boundary 
conditions possesses 2(L — r) largest (in magnitude) eigenvalues which are asymptotically 
degenerate in the thermodynamic limit, corresponding to 2(L — r) coexisting phases. These 
largest eigenvalues are separated from the other ones by a gap which remains finite in the 
thermodynamic limit. The corresponding 2(L — r) (quasi-) ground states can be constructed 
by means of Bethe Ansatz. In this framework, the spontaneous staggered polarizations can be 
computed, as in the case of the .F-model pQ, from the normalized matrix elements, between 
two of the corresponding Bethe eigenstates, of the a^ Pauli spin operator acting on the spin 
(or variable e m = ±1) on a given bond m of the lattice. 

The computation of spontaneous staggered polarizations in the CSOS model has already 
been considered in [2] in the case 77 = 1/L. The derivation of [2] uses the representation 
of Bethe vectors in the framework of coordinate Bethe Ansatz. However, since no compact 
and convenient representation was known at that time for the scalar products (and form 
factors) between such Bethe eigenstates, the derivation of [2j actually relies on two conjectural 
mathematical identities for elliptic functions. The use of the determinant representations 
obtained in [9] from algebraic Bethe Ansatz enables us instead to easily compute (for any 
rational value of the parameter 77 of the model) the thermodynamic limit of the form factors 
contributing to spontaneous staggered polarizations. By this method, we are also able to 
control the finite-size corrections to the result in the same way as in [TJ. 

This paper is organized as follows. In Section [21 we briefly recall the ABA solution of the 
CSOS model and the result of |9J concerning the finite-size determinant representation for 
the normalized matrix element of the local operator a^ between two Bethe eigenstates. In 
Section O we characterize the 2(L — r) (quasi-)ground states from the study of their Bethe 
equations, and discuss the process of taking the thermodynamic limit and of controlling the 
finite-size corrections. In Section [H we apply this process to the form-factor representation 
of Section [2] in the case where the two Bethe states belong to the previously described set 
of ground states. We obtain a Fredholm determinant representation that can be explicitly 
computed, leading, in Section [5l to explicit expressions for the spontaneous polarizations of 
the model. 

2 The form factor in the ABA framework 

Let us consider a two-dimensional square lattice of size N x N (JV even), with periodic bound- 
ary conditions. The statistical weights H / ( o , s f +{E ? ), e^, e',e 7 ,e' € {+1, —11, corresponding 



to the six allowed configurations around a face of the SOS model, can be seen as the six 



non-zero elements R(ui — £j; s) V ] of the following i?-matrix: 



R(ui -£j-,s) 





•ajJ-/ 


Wi 












/! 
















°\ 





b(iti 


~ti 


s) 


c(Ui 


-e, 


s) 








c(v,i - 


-ir, 


s) 


b(v,i - 


-fc; 


s) 





\o 
















V 



GEnd(C 2 ®C 2 ). (2.1) 



Here Uj (respectively £j) is an inhomogeneity parameter attached to the column i (resp. row 
j) of cells of the lattice. The functions b(u; s) and c(u; s) are given as 

b(u-,s)= [ * + 1] M c(u;s)= [s -\ u][ ] ] with [u]=ei( V u;r), (2.2) 

[s\ [u + 1J [sj [u + 1\ 

where 9\ denotes the usual theta function (jA.lj) with quasi-periods 1 and r (9r > 0). The 
height s in (|2.1|) is called dynamical parameter. The i?-matrix (|2.ip with dynamical parameter 
s satisfies the dynamical Yang-Baxter equation [61 [3], and the corresponding SOS model is 
also sometimes called dynamical six-vertex model. In the cyclic case that we consider here, 
the parameter t\ of the model is chosen to be rational: rj = r/L, r and L being relatively prime 
integers. Hence, the statistical weights elements of the i?-matrix are L-periodic functions of 
s. 

Remark 2.1. Our parameter rjso, where sq is a global shift of the dynamical parameter (such 
that s — so is an integer) introduced so as to avoid the singularities in (|2.2p , is related to the 
phase angle fi = wo/ir of the physical model introduced in [121 E] by a shift of r/2 (later on, 
we will for simplicity take so = •£•, i.e. /i = 0, in our result). Also, the statistical weights of 
[12} [TT] correspond to a diagonal dynamical gauge transformation of the i?-matrix (|2.ip which 
leaves the local height probabilities invariant. 

Eigenstates of the dynamical transfer matrix constructed from the .R-matrix (|2.ip can 
be obtained by means of algebraic Bethe Ansatz (see [US]). In this framework, the space of 
states of the model corresponds to the space of functions Fun(%[0]) of the dynamical parameter 
s € so + TLjLTL with values in the zero-weight space "H[0] = { |v) G % \ (J2j=i a j) l v ) = 0}' 
with % ~ (C 2 )® N (see [9] for more details). More precisely, from a product of i?-matrices 
along a column of elementary cells of the lattice, one can construct the monodromy matrix, 
which is a 2 x 2 matrix of operators 

?(„)=(*■> |'">) ,2.3) 

\C(u) D(u)J 

with A(u),B(u),C(u),D(u) € End(FunH). Common eigenstates to transfer matrices t(u) = 
A{u) + D(u) can be constructed in Fun("H[0]) as 

n [11 
| {v},u ) : s m- oo s TT ^^ B(yi; s)B(v 2 ; s - 1) . . . B(v n ; s - n + 1)| 0). (2.4) 

f=i t s " 3] 

Here n is such that A^ = 2n + NL for some integer K, u is such that (— l) rn uj = 1, the 
reference state | 0) S % is chosen as | 0) = <8>£ 1 ( ), and {v} = {v\, . . . ,v n } (with rjVi ^ i]Vj 
mod Z + rZ) is a solution of the system of Bethe equations 






j = l,...n, (2.5) 
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Similarly, eigenstates of the transfer matrix in the dual space are constructed as 

if 

[1] 



( {«}, w I : s M- ( |CK; s - n) . . . C( W2 ; s - 2)C(v 1 ; s - 1) uT s \{ ±^I, (2.7) 

with ( | = | )t and {v} solution of (f23j) . Then 

?(u) \{v},u) = t(u;{v},uj) \{v},uj), ({v},u\t{u) = t(u;{v},uj) ({v},uj |, (2.8) 

with 

r(« ; w, W ) = u, »(«) n [ t m+ i 1] + (-i)***- 1 « n [ w ] - ( 2 - g ) 

AA [ v , _ u j XX [u-Vl\ 

From now on, we restrict our study to Bethe states in the sector n = N/2, which is 
enough for our purpose. For two such Bethe eigenstates, associated to two different solutions 
{u},to u and {v},co v of the system of Bethe equations (|2.5p . we want to compute the following 
renormalized form factor: 

s m ({u}, uj u ; {vj, w„) = — yj (2.10) 

{({u},cj u \ {u},u u )({v},uj v \ {v},U} v }) 

{{u},uj u \a^\{v},uj v ) f ({v},uj v | {v},uj v )\ 1/2 



( {v}, bJ v | {v}, U) v ) V ( {">}, ^u I {u}, U} u 

This quantity can be rewritten in terms of a ratio of determinants by means of the represen- 
tations obtained in our previous paper [9]. We recall that the matrix element of the operator 
cr^j between two different Bethe eigenstates {{u},lo u | ()2.7[> and | {v},co v ) (|2.4|> in the sector 
N = 2n is given as 

II,.\ ,, \„* \U,\ ,, \ f TT r (&'>M» w «) ] Jl V- ..-*..» [7 + «] 

( {«}, W u | (7 m | {V}, W„ ) =" 






'' l ' /l " / ' ■det[fi 7 ({n},o; n ;{«},a;„)-2P 7 ({n},a; n ;{w},w t ,|C m )], (2.12) 



Wk<i\ u k ~ ui][vi - v k ] n 



with 7 = Y^ij=x( v 3 ~ u j)- Here f2 7 is the matrix appearing in the determinant representation 
for the scalar product (see Theorem 3.1 of [9]), with matrix elements 

rn ti i r i Yl X /[ u j -u *+7] u v [uj - v k + j + 1]\ -A- 

[^(M,^;W,a; v )]. fe = H |^— ^ [u ._ Vk + 1] ^(v k )^[u t - v k + 1] 

I {[Uj-Vk+J\ U u [Uj - V k + J - l] \ _ 2 ,/ A TTr rl , orQ x 

+ T"-r< — 7 i 7 7Tr w « o(wfc)l Ipt - Vfc - A > 2 - 13 

7 [ «j - ^fc w,, «j - w* - 1 J f_\ 



whereas V^ is a rank 1 matrix denned as 



r-n rt i r i imI X / K ~ &»» + H w„ [uj - £ m + 7 + 1] 



jk [7] I [«J - Cm] W M [«j - Cm + 1] 

We also recall the determinant representation for the 'square of the norm' of a Bethe eigenstate: 

< WlM . 1 WlM . ) = n;., -hw <m^H -» + '] det [#(M)] , (2 , 15) 

(-[oj') n ll^feK- -^J 

with 

A ^ [«j -u t - 1]' [Uj -u t + I}' 



[•«»})]„ = «J log- 2 w + s -"'-*',' - ■"'-"' + !',' ) 

Jjfe [ d ^\K-«t-l] [«i-t*t + l]/J 

[«j ~U k - 1]' [lij - Ufc + 1]' 



[Uj - u k - 1] [«j - u k + 1] 



(2.16) 



Factorizing the quantity a(vj) YYt=i[ v t — Vj + 1] = ~^ v 2 d(vj) TYt=i[ v t ~ v j ~ 1] ou t °f each 
column of the determinant in (|'2.12p . we can rewrite the first ratio in (|2.1ip as: 

{{u},uj u \a^\{v},uj v ) )n r(4;{«},w u )| J 1 ^ /^V^ + s]! 



n T(£k){v>},Uu) I I 1 V^ 



{v},u v I {v},u> v ) I J-l r(C fc ;{v},a;^) I L ^ /r _ W 



fe=i 



■ n1 An TT d K) TT N ~ ^] det ra [ff({^}, {v}) - 2Q({u}, {v})] 

' ' l \ d(v k ) 11 K " «l] det n [cD({t,})] ' l • J 



where 7 = X)t=i( u * ~ u t), ^ is given by f)2.16fl . and 



rHY/,,1 /„\M - 1 ( ^ ~ gj ± 7 ^ ^ [^ ~ V A ± 7 ± 11 1 TT [^ ~ ^ + X J 

L*U«M«»Jy- [7] j [u| _ Wi] ^ [„,_„, + !] jl = l[ Ut - Vi + i] 

1 / [u, - uj + 7] u u [m - Vj + 7 - 1] | cj -A- [u« - uj - 1] 

[7] J [Ui - Vj] U V [Ui - Vj - 1] J w2 11 [^ - „j - 1] ' 

[n/ ; w,n, 1 f [Ui ~ Cm + 7] Uy [Uj ~ Cm + 7 + 1] ) A [^ ~ Cm + 1] ,„ n Q x 

In its turn, from (|2.15p . the second ratio is given by 

( {v},u v | {v},u v ) _ -pr a(v k )d(v k ) jT \ v i - v k + l] tt K' - u k\ det n [$(>})] _ 2Q . 
({u},u; u | {u},u; u ) -J Jra^*)^*) ,_Auj - u k + 1] *J:[vj - v k ] det n [$({u})]' 

These representations are the starting point for our study of the spontaneous staggered polar- 
izations of the model, which can be obtained from the thermodynamic limit of the quantity 
(|2.10p in the case when | {u},uj u } and | {v},uj v ) correspond to two different (quasi-)ground 
states of the homogeneous model. 



Remark 2.2. The case of the mean value of the operator a^ in the same Bethe state should 
be treated separately since the proper limit has to be taken into the determinant of (|2.12|) in 
the same way as in (j2.15j) - f|2.16[) . We will in that case use the formula 

„*„,-, iX , det n [$(M) + 2Q(+)(M)] 

s m {{u},uj u ;{u},u u ) = , (2.21) 

det n l<P({u})\ 



with $ given by (|2.16p and 

\0 (+) (lu\)] - ^ Uj - ^ - ^ Uj ~ ^ m + ^ (2 221 



One should also pay special care to possible other cases for which 7 = or when some of the 
Vj coincide with some the Uj. In particular, when L is even, there is a little subtlety that 
what not mentioned in our previous article [PJ: in that case, with each solution {u},u 2 of 
the Bethe equations (|2.5p . one can associate two different Bethe eigenstates |{u},w) (with 
u = e to(rn+2t)/L for some integer £) an(1 | { u } ; _ w ) (_ w = j^rn+wyL w ith e = e-L/2) with 

opposite eigenvalue. The corresponding normalized form factor can be explicitly represented 

as 

*(i x s x ^ (~ 1 ) m " 1 V- -*™ det w [^-)(M) + 2Q(-)( 8 )] 

s£s +Z/LZ n L Vl J; J 



with $ given by <|2.16j> and 



[*'->(W)]^<4<W + |(|^I 



ut+ir 



[Uj - lij + 1] 



+ , [uj -nfc- 1] / _ K-ufc+ir , , 2 24 ^ 

[Uj - Ufc - 1] [Uj - Ufc + 1] 



fsl' 



[Q H WU = 2^- (2-25) 

3 The degenerate ground states in the thermodynamic limit 

According to |12|. II lj . the ground state of the CSOS model is degenerate at the thermodynamic 
limit: there are 2{L — r) (quasi-) ground states in the sector n = N/2, which were identified 
in |12l [TT] in the low-temperature limit. 

To characterize these states, it is convenient, by means of Jacobi's imaginary transforma- 
tion 

[u] = -i(-zr)- 1 /2 e ^w« 2 01 (^ ;f ) ) f = --, fj = --, (3.1) 



to rewrite the Bethe equations (J2.5J) in terms of theta functions with imaginary quasi-period 
f . In the sector n = N/2 and at the homogeneous limit £& = 1/2, k = 1, . . . , N, it gives 



^(^-f?/2)l = l^-^-f7) X ' ^ 

where we have set Zj = fjvj, j = 1, . . . , n. Here and in the following, unless explicitly specified, 
the considered theta functions are of imaginary quasi-period f, i.e. 9\(z) = 8i(z;f). 
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These Bethe equations can be rewritten in the logarithmic form as 
n -. n 

N Po (z j )-Y,^j-zi) = 27r(n j - 1 ^- + (3 + 2r ] Y,zi), j = l,...,n, (3.3) 

1=1 1=1 

where rij are integers, po and •& are the bare momentum and bare phase 

Po{z) = llog 9 1 (f l /2- z y m = llog eJ?Tz)' ( ) 

and w = e i7r/3 . The degenerate ground states identified in [TT] correspond to real solutions of 
(|3.3p such that %+i — rij = 1, j = 1, . . . , n — 1. In the thermodynamic limit N — )• oo (with 
n = N/2), the distribution of the Bethe roots corresponding to such states tends to a positive 
density p(z) on the interval [—1/2, 1/2], solution of the following integral equation: 

p(z)+ f K(z-w)p(w)dw = ^-, (3.5) 

7-1/2 2tt 

with 

K[Z) -2n (Z) -2 7 r\0 1 (z + rj) O^z-i})}' ( j 

The solution of the integral equation (|3,5p can easily be computed by means of the Fourier 
transform. In the domain < r/ < 1/2, i.e. < — ifj < — |r, we have 

7 _ I ,1/2 ^ .-a**™ ^ form = 0, 



p' m = p' (z)e- 2mm *dz=t . Hl .jM r H + , . (3-8) 

7-1/2 [2vrg 2 ^_^ m| — otherwise, 

/■V2 _ . fl for m = 0, 

fc m = / K{z) e~ 2mmz dz = Hm| 1 _piH g -2|H f , . (3.9) 

7-1/2 [<? i_p| m | — otherwise, 

such that 

ts_ v^ e + tt I 1 -g J V3(z,V) mn \ 

P{) ^ 2cosh(i7rmr?) 2 H (1 + q m f 0Az;fj) ' l J 

m=— oo v " m=l 

in which we have set q = e 2mfl , p = e 27nr . 

We now want to study more precisely how the Bethe roots of one of these ground states 
behave with respect to finite size corrections. Let us introduce, for a given (quasi-)ground 
state parametrized by the set of roots {xj}i<j< n solution to f|3.3[) with uj x = e l7rl3x and a given 
shift of integers k x = rij — j (which does not depend on j), the following counting function: 

1 1 n 1 / -I- 1 n \ 

Uz) = -p (z) -^J2^ z - x ^ + -[ ! V~ - k* - A* - 2tj J>l ) ■ (3.11) 

71 n i=i n V 1=1 J 

This function is such that £ x ( x j) = i/ n ) j = 1> • • • > n - Moreover, since its derivative, 

1 2 " 

£(*) = -p£,(s) - - ^ K(z - Xl ), (3.12) 

j=i 

tends to 2p(z) (which is positive) in the thermodynamic limit, £ x is an increasing, and hence 
invertible function, at least for N large enough. Hence one can show the following result: 



Proposition 3.1. Let f be a C°° 1-periodic function on R. Then, the sum of all the values 
f(xj), where the set of spectral parameters {xj}i<j< n parametrizes one of the degenerate 
ground states solution to (j3.3|) . can be replaced by an integral in the thermodynamic limit 
according to the following rule: 

i n r 1 / 2 

J7 E K x j) = / /(*) P^ dz + 0{N~°°). (3.13) 

Similarly, if g is a C°° function such that g' is 1-periodic, then 

• 1/2 



1 P _1_ / A 

^E^i)= / 9{z)p{z)dz+ C j-Y,x ] +0{N- 

.7=1 • 7_1/2 .7=1 



(3.14) 



-1/2 



^ere c ff = J_i/ 2 </(*) tte = 5 (l/2) - ^(-1/2). 



Proof -- The proof of ()3.13[) is similar to the proof of Proposition 3.1 of [7\. For completeness, 
we recall its main arguments. It relies on the fact that one can easily prove, using the Taylor 
expansion of the 1-periodic function /, the analog of (|3.13p in the case of homogeneously 
distributed variables: 

-E/(-)= / mdz + 0(N-°°). (3.15) 

One should then notice that £ x is a C°° function of real variables such that ^(z+l) = £ x (z) + l, 
so that the function / o ^~ l is also 1-periodic. One can therefore apply (|3.15j) to / o ^~ l and 
perform a change of variables in the integral to express the sum over Xj as 

-E/(^) = -E/fe'^'/n)) = / f(z)&(z)dz + 0(N-°° ). (3.16) 

n i=i n ;=i - 7 - 1 / 2 

Finally, applying (|3.16|) to the r.h.s. of (|3.12|) . we obtain that, up to corrections of order 
0(N~°°), the function £ x /2 satisfies the same integral equation (|3.5p as p. By unicity of the 
solution, we have 

Z(z) = 2p(z) + 0(N-°°), (3.17) 

which ends the proof of ()3.13p . 

The identity <|3.14|) is then a direct corollary of (|3.13p : if g'(x) is 1-periodic, then g(x) — c g x 
is also 1-periodic, and one can apply f|3. 13|) to get 

^E^)-iE^'= / 9{z)p{z)dz-c g zp(z)dz + 0(M-°°), (3.18) 

j—l j—i J— 1/2 J— 1/2 

the last integral being zero by symmetry. □ 

Remark 3.1. Proposition 13.11 can be used to obtain a sum rule for the corresponding ground 
state roots. Summing all logarithmic Bethe equations (j3.3|) for j = 1, . . . , n and using the fact 
that ■& is an odd function, we get 

1 n n 

77 i=i j=i 



Using then (|3"7T4"1) applied to the l.h.s of (|3"7L9i 

n n 

£>(*,-) = 2vr J>, + 0(N-™), (3.20) 

we obtain that 

X> = ^| + O(iV-) = L ^™_ + ** + Q^- 00 ), (3-21) 

in which we have set ^ = rn + *» . 

Remark 3.2. The counting function (|3.1ip can be evaluated in the thermodynamic limit as 

l(z) = 2 f p(w) dw + ^±l-?-J2 x i + 0(N~°°). (3.22) 

This follows from (|3.17p and from the value ^(0) which can be evaluated in the thermody- 
namic limit by means of (|3.2ip as well as (|3.14p applied to the odd function §: 

1(0) = - ;4l>(-^) + U 1 ^ 1 - k. - A, - *,Y,x 3 ) 

1=1 V 3=1 J 

One can now use these results to compute more precisely the infinitesimal difference of 
roots X7-4.1 — Xj. The latter is given as the value S x (xj) of the infinitesimal shift function 5 X 
defined as 

5 x (z) = £- 1 (Uz) + l)-z- (3-24) 

Rewriting the equation ^ x {Sx(z) + z) = £, x (z) + ^ using the representation (|3.22p of £ x , one 
gets 

z+5 x (z) -1 

p(w)dw = — + 0(N-°°). (3.25) 

Expanding p in Taylor series, one therefore obtains a relation which enables one in principle 
to compute 5 x (z) at all order in N: 

00 

N H T\ p{k ~ l){z) ^ x{z) ^ = 1 + °( N ~ CO )- (3- 2 6) 

k=l 

Let us now consider two different ground states for the system of Bethe equations (|3.3p . 
parameterized by a solution {xj}j=i,...,m ui x = e l7T ^ x , and a shift of integers k x (respectively 
by {yj}j=i,...,n, u y = e l7T P y , and a shift of integers k y ). We want to evaluate, at large N, the 
infinitesimal difference of roots Xj — yj. To this aim, we define the infinitesimal shift function 

&,„(*) =ex X (&(*))"*, (3-27) 



which is such that 6 x ,y{yj) = %j — Vj- Using again the representation (|3.22p for S, x and £ y , one 
obtains for 5 x ,y an equation analog to (|3.25j) : 



/ " p(w)dw = -Y,(x3-yj) + 0(N-°°), (3.28) 



which leads to 

oo 1 

N H u p{k ~ l) ^ [£>*(*)]* = 5>' " Vl) + 0(N~°°). (3.29) 



kV 

k=X 



Comparing this equation to (|3.26[) , we obtain that the two infinitesimal shift functions 5 X and 
6 x>y , seen as functionals of the density p, are related by 



$x,y [p( z )] 



p{*) 



Y,ii x i-yi) 



0{N-°°). (3.30) 



This characterizes, at all orders in N, the infinitesimal difference of roots Xj — yj. We recall 
that, from (|3.21|) . the sum Yli( x i ~ Vl) ls itself given by 

I> " Vl) = K *-^ + _ P *~ p y + 0(N~°°) (3.31) 

(3.32) 
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*-y ~r P X 
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+ 0(N- 


-oo\ 
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~ y) ' ^K^x 


~Zy) , 


0(N 


at 


2{L- 


r) 
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Vl) — 


e i-K(k x -k y 


) w x 

UJy 


+ 0(N- 


-oo\ 



which also means that 

e 2m(i-T))J2i(*i-vi) — g^i.Ktc-Kyj ~z± _|_ 0(N'°°). (3.33) 

LOy 

In particular, the previous study enables us to identify (and count) the degenerate ground 
states associated with a given set of Bethe equations. We see from (|3.30p that two solutions 
{x} and {y} are different if and only if the total shift Yli( x i ~ Vl) ls n °t an integer. Hence, in 
this setting, the different degenerate ground states are completely determined by two quantum 
numbers k € %jTL and £ G Z/(L — r)Z. If L is odd, this gives 2{L — r) different values (modulo 
1) of (|3.32|) . corresponding to 2(L — r) different Bethe eigenstates. For even L, we only get 
{L — r) different values of (|3.32[) . each of them being associated to two opposite values of oj 
according to the parity of k (see (|3.33p ). 

4 The form factor in the thermodynamic limit 

We now study the thermodynamic limit of the renormalized form factor (|2.10|) in the case 
where | {u},uj u ) and | {v},uj v ) are two ground states of the homogeneous model. 

Let us set Xj = fjuj, yj = fjvj, j = 1, ... ,n, and u x = u u , uj y = oj v . We have seen in 
Section [3] that the Bethe roots Xj (respectively yj) for one of the ground states are completely 
determined by the data of two quantum numbers k x ,£ x (respectively k y ,£ y ). From now on, 
we simply denote \k x ,£ x ) = \{u},u u ) (respectively \k y ,£ y ) = \{v},u v )) the corresponding 

Bethe eigenstate, and | ^g x ' x ) (respectively | tfj g y ' y )) the corresponding state renormalized 
to unity, i.e. 

\^.(k x i x )\_ I kg, 4) |,7.(ky.M^_ \ k y^y) ^^ 

\\k X ,l X \ k X ,l X }) ' { ( Y.y , ly | Ky , ly ) ) ' 
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In order to study how the determinant representation for s^ behaves in the thermodynamic 
limit, it is convenient to rewrite it in terms of theta functions with quasi-period f by means 
of Jacobi's imaginary transformation (|3.ip . Let us first consider the quantity (|2.17|) in the 
case 7^0. We obtain 

{KJx\<7m\K>ty) _ ^„-27mn -i£" = ibo(>i)-po(y0 



ky,iy | ky,ly ( 



-e — -"e 






tt gi(a/fc - W) det n [H({x}, {y}) - 2Q({g}, {y})] 
lie^xk-xi) det n [$({y»] 



with 7 = YJj=i{yj - x i) and 



r^n 1 f gi (gj - y fc + 7) u; y ^^ 0i ( gj - - y k + 7 + 77) \ -A- ft (x; - y k + 77) 

L J ^ #1 (7) I #1 (xj -Vk) u x d x (xj -y k + fj) J Al 0! (y, -y k + fj) 



1 f gjXgj - j/fc + 7) _ ^£ -2irirry 0l( x j -Vk+l-rj) 



n\ 



e 



(7) \ 0i(xj-y k ) uj y 9i(xj - y k - fj) 



x 



\ m_ / J--1 Hi 1 1/7 — 1/7. — ri I 



^/ ,J;Oi(yi-yk-fj) 



1=1 



[gl 1 [ g l(Xj - -| +7 ) a, ^ AxQej +1+7) lfr gifo + |) f44 , 

[*] ifc = -2nif}NS jk t?^ - ^Y,K( yj -yi)\- 2irifjK(yj - y k ) + 47vifjr). (4.5) 



The determinant of the matrix <J> (|4.5|> is already in a convenient form for taking the 
thermodynamic limit. Using (|3.13p . (|3.5p . we get 

[5], t . -2«*) {* + i^g^ - i^j + <**-)} , (4.6) 

so that the corresponding determinant can be written in terms of a Fredholm determinant in 
the thermodynamic limit: 

n 

det[$({y})] =(-2vrir/iV)™J]p(y,){det[l + ^-Fo]+0(A r " 00 )}. (4.7) 



1=1 



Here .ff and Vo are integral operators acting on the interval [— ^, o]j with respective kernels 
K(y — z) given by (|3.7p . and Vo(y — z) = 1r\. 

The determinant appearing in the numerator of (|4,2|) can be transformed, similarly as 
what was done in [7], in a more suitable form for the thermodynamic limit. Using the results 
of Appendix [Bl we have 



i(Vj ~Vk) | 

-L-{det[(H-Q)({x},{y})] - det[(W+ C)({x},{y})] }, (4. 



(- 1 )" II a (J f d „ et [^(W. M) - 2 «(^}' {V})) 

j<k 1[ j k) 
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with 

l n \jk °jfc n n „/ _ n S 11 &(.„_,„ ,fi\ 1,,>J LLf), 



]ir=i i (Vj - x i)\fJ[6i {Vl -yk + v) ^ u x J jLi 9i (yi -yk-fj) 

1 u y jbtirn f O'liVj ~ Vk + V) _ O'liVj - yk - fj) 



-I ?-e 

9[ (0) w x [ 6i (yj -yk + v) Qi (Vj -Vk-fj)) ' 



5 jk Nct>-\{x},{y}) U + (yk\{x},{y}) ~ (^) 2 e 4 ^^-(y fc |{x},M) 



27ri w. 



".y 



?i(0)w a 



e 2 ^K( yj -y k ), (4.9) 



and 



rpl = ^. e 2 ^7 _ rr e i( x i + V °i(yi f) = ^2™^ _ -iEr=ifeo(*j)-po(M)] f4 101 

In (|4.9p . the factor cpj({x}, {y}) and the functions 4>±(y\{x}, {y}) are respectively defined as 

The behavior of the matrix elements of Q can straightforwardly be evaluated in the ther- 
modynamic limit. Prom (|3.20p and (|3.33p . we obtain 

[Q] jk = [(-l) k -l]e 2 ^ + 0(N-°°), (4.13) 

where we have set k = k y — k x . It follows in particular that the form factor ()2.10p vanishes in 
the thermodynamic limit when k = 0. 

The behavior of the functions <p±{y\{x}, {y}) (|4.12p . conveniently rewritten in the form 

0±(y\{x},{y}) = exp< ^\tog0i(x l -y±fj)-]osO 1 (yi-y±fj)) >, (4.14) 



i=i 



can be evaluated by means of Proposition [370 Indeed, the functions g^ : z i-^- log 9i(z — y±fj) 
being C°° on R with a 1-periodic derivative, one can apply (|3.14p to each of the sums in (|4.14p . 
with corresponding constants given by 



c g± = J n£-JL±l»dz = ^. (4.15) 



'- 1 ' 2 6' 1 (z-y±fj) t 
.1/2 6i(z-y±TJ)' 

One obtains 

</>±(y|{*}, {y}) = e ±l ^ + 0(N-°°). (4.16) 

So as to evaluate the factor (j) (|4.1ip . let us define the function 

12 



which is C°° and 1-periodic on R. This function is such that 

#Vilto,{i/}) = -^|^^(W,{y})- (4.18) 

Using the fact that the sum of the residues of the elliptic function 
0\{z-y) f-J^ e^z - yi ± 7]) 



vanishes in an elementary cell, we obtain the following identity for the quantity 



"3- 



jf £ ^(y-yjllj® *'«*>' {y}) = ^ ^ (yl{XJ ' {y}) - (420) 

Rewriting the summand in terms of arguments of 1-periodic C°° functions evaluated at the 
ground state roots yj by means of (|4.18p . and applying Proposition 13.11 we obtain that, in 
the thermodynamic limit, the function (ft (|4.17p satisfies the integral equation 

jm. \ V ' "f-'+jy m i; = #« e^) + 0(N-~). (4.21) 

2iir J_ l/2 6i(y-z±r]) 

The latter can easily be solved by Fourier transform and we get 

-4^Wi) = 4>i pfa) + 0(Ar-°°) = sin(7T7)^^ + 0(N~°°). (4.22) 

Gathering all these results, and using also (|3.33p . we can therefore rewrite the elements of the 
matrix %({:c}, {y}) (|4.9p in the thermodynamic limit as 

m* = -^§> ^ {** + ^^r + ° (iNroo) } • (423) 

Hence, the quantity (|4.2|) admits the following representation in terms of Fredholm deter- 
minants: 

(k x ,e x \a^\ k y ,£ y ) _ ( nk ( m _i)J 1 ^ (^y_ e 2irirfyY 9i{f\s + l) 



( _ 1)k(m -D ii ^ 



KAWv) " ' l L seso ^/LZ V ^ ^i(^)^i(7) 



L)K _., , det[l + (-l) k K + ^^F]-det[l + (-l)^-^^i>] 



^x 



u y J det[l + K-V ] 

+ 0(N-°°). (4.24) 

We recall that K has kernel K (|3.7p . whereas Vo and V have respective constant kernels 2r] 

and^i(O). 

In its turn, the second ratio (|2.2Up can be evaluated using the Bethe equations for {u} and 
{v}, the estimations (|4.16|) and (|4.22p . as well as the Fredholm determinant representation 
([477]) . It gives 

( ky, A/ i k ^ *v ) ._ (<\ U + (N-°°), (4.25) 



{k-xj^x \ k. x ,l x } \U 
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so that the conveniently renormalized form factor (j2.10|) is simply given by 



(4 kM \- Z m \4 yM ) = (-l) k(m ' 1) \j. £ 



L 



L (r?s + 7) 



S£SQ + 



w* / 0i (77s) 0i (7) 



det[l + (-l) k if + ^Y' V] -det[l + (-l) k ET 
det[l + ^-T> ] 



i-(-i) k 



V] 



+ 0(N-°°). (4.26) 



The above Fredholm determinants can be computed by means of Fourier transform. In- 
deed, as the kernel of the integral operator K depends only on the difference of two variables, 
its eigenvalues correspond to the Fourier coefficients k m (|3.9p of the function K (|3.7p . We 
obtain 



det 



[1 + K-V ]=(2-2 V )Y[ (l + «T— 



rra=l 
hoo 



I — ff n Q- 2m \ 2 



P" 



2(1 - r,) J] 



m=l 



(l + q m ) 2 (l-p m q- m ) 
(l-p m ) 2 



(4.27) 



whereas 



det 



l + (-l) k K± 



1 



•l) k 



V 



1 — (— ~\) w 
! +(_l)k±i ^ ^^i(0) 



2tt 



n 

m=l 



{l + {-l) k q m ) 2 (l-(-l) k p 
(l-p m ) 2 



k~mp;—m\2 



(4.28) 



Finally, the renormalized form factor (|2.10p is given as the following infinite product 

1 l 1 \ k +°° ft ~m\2(-\ 1 ~m~—m\2 

(^ M 1 ^ 1 4 vM ) = ( } (-1)"" 1 TT g 9 +p g 



7r(L — r) 



• (1 + q m ) 2 (l -p m g- m ) 2 
y eMa ?-i). ^ + 7) fl(0) 1 Q 



sgso+Z/ 



The previous expression is a priori only valid in the case 7/O. One can perform a similar 
study in the two particular cases mentioned in Remark 12 .2\ i.e. for the mean value (|2.21|) and 
the form factor (|2.23p in the case L even. One obtains that the former vanishes, whereas the 
latter can be written as 



<v> 






) = (-1)— 1 11 |_ 



(l-^ m ) 2 (l+P 



mz—m\2 



IC 



-17TS0 



• (1 + q m ) 2 (l -pmg-my 

B-vfj^t-J+or"). (4-30) 



m 
L-l 



7r(L — r) 
Note that we can rewrite a general representation valid for all cases by conveniently reg- 
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ularizing the representation (J4.29P as 

( ^m ]o m ,/,«)_ i-(-i) k ( ir -i ff ( 1 -n 2 (i+P m n 2 



x lim 



y e i,^-i +27?Q)s ft(fr + 7 + a) e[(0) 1 Q 



o I 7r(L - r) 

We recall that, from (|3.32|) . the difference of roots 7 can be expressed in terms of the quantum 
numbers k = k y — k x , £ = £ y — £ x as 

TV _i_ op 

7 = ^^ + 0(iV- 00 ), (4.32) 



so that 



{ ^m 1 a > 1 j»*m ) - i-(-i) k f ir i fr a-n 2 (i+p m -—' « 2 



f (1 + <f ?1 ) 2 (l -p m £- m ) 2 

-(^+27^ + ^1 






seso- 

(4.33) 

Note that this quantity depends only on the differences k and £ of the quantum numbers: 

(^ A) I < I 4 WM ) = *4(M) + 0(N~°°). (4.34) 

5 Spontaneous staggered polarizations of the CSOS model 

Let us consider the 2(L — r)-dimensional subspace Fun(% 5 [0]) of the space of states Fun(%[0]) 
generated by all the ground states, i.e. by the Bethe eigenstates associated to the 2(L — r) 
largest (in magnitude) eigenvalues of the transfer matrix in the thermodynamic limit. A basis 
of this subspace is given by the normalized Bethe eigenstates \ip g °" a ), with k a G Z/2Z and 
£ a G Z/(L — r)Z. This basis is not polarized, since the mean value of a^ in one of the Bethe 
ground states \ip g °" ) vanishes in the thermodynamic limit. A polarized basis of Fun(% 9 [0]) 
is instead given by 

I 4' t} ) = ,J l E ^(-l^'e-***^^ I ^M ), (5.1) 



# " r) 



k Q =0 f a =0 



with e G {0, 1} and t G {0, 1 . . . ,L — r — 1}. Spontaneous staggered polarizations of the CSOS 
model are given as the mean values of the a^ operator in the polarized states (|5.ip . 

Remark 5.1. It can be shown, for instance by considering the combinatorial formula for the 
Bethe states (see Theorem 5 of [1]), that the combination (|5.ip tends, in the low-temperature 
limit p, q — > 0, to one of the flat ground state configurations of the type (a, a + 1, a, a + 1, . . .) 
or (a + 1, a, a + 1, a, . . .) identified in [12] • More precisely, in the case so = £ = — <4, the state 

I 0p } with t = a — [r]a\ tends to the flat ground state configuration (a, a + 1, a, a + 1, . . .) if 
e — L^aJ is even, and to (a + 1, a, a + 1, a, . . .) if e — |_??aj is odd (here |_^J denotes the integer 
part of x) . 
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In fact, the matrix elements of the a^ operator in the states ()5.ip are given as 

(4 ei,tl) l<l4 £2 ' t2) ) = ^7 I L 7y £ £ (-i) kl61+k2£2 

V ' ki,k 2 =0 £i,£ 2 =0 



■[rh±^ { t 1+so )-r^±^ {t2+so) ] { ^ (kli/l) ! ^ ! ^ {k2i , 2) ) 



xe L _ NY , 9 

1 1 L-r-l 

1 \^ V^ C_i^k 1 (ei-e 2 )+k£ 2 

2(L-r) ^ £-" y ' 

v ; ki,k=0 €i,€=0 



27T 

x e 



[rki+2^1/ > rk+2^/4. , „ \1 

L^^- (tl ~ t2) -^^ (t2+so) J s z m (k,£) + 0(N-°°), (5.2) 



in which we have set £ = £2 — ^1 and k = l<2 — kj.. We see from the sum over t\ that this 
quantity is non-zero only if ti = t2, and then from the sum over ki that one should also have 
ei = 62- Hence, up to exponentially small corrections in the size of the system, the local 
operator o z m is diagonal in the basis ()5.ip . 



( 4 £l,tl) I < I 4 2,t2) ) = Wa ^,t 2 ( 4 lM) I o z m I 4 ei ' tl} ) + 0(JV-°°), (5.3) 

and the corresponding mean values are given as 

•/ -i\m-l+e +°° (1 ~m\2(i , ~m~-m\2 

Wfl I ml^ 9 / Tr(L-r) J- 1 (l + g m ) 2 (l-p m r m ) 2 

m=l 

The expression (j5.4|) coincides with the one obtained by Date et al. in [2] in the case 
r = 1, L odd and so = j- = — n^- Note however that the result (|5.4|) can be presented in a 
much simpler form by means of the identities of Appendix [A] for the sums of theta functions. 
Indeed, using respectively (|A.6|) to compute the sum over £, and then (|A.7|) to compute the 
sum over s, we obtain a simple combination of theta functions in which the a — > limit can 
be straightforwardly taken, so that 

TT 6*2(0; r/rj 6*2 (r/t; (1 - r/)rj 

^ (l-g™) 2 ^-^-™-*) y? (l-p m q- m+t ) 
11 (1 + guy (i + pm^-m-t) 11 (1+^-m+t) 

+ 0(iV-°°), (5.6) 



•1 



.m— 1+e 



in which we have set sq = —j=- 

One recovers the low-temperature completely ordered limit by taking |f | — > oo, i.e. p,q — > 
in this expression. If on the contrary one wants to study the critical limit \t\ — > oo, it is 
better to re-express the result (|5.5p by means of Jacobi's imaginary transformation as 



( <T } I < 1 ^ ) = (-D m+£ | ^o-^l^-M + ° (Ar_00) - (5 - 7) 
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A Theta functions and useful identities 

In this paper, 6\(z; r) denotes the usual theta function with quasi-periods 1 and r, 



oo 



k=— oo 

which satisfies 



1 (z;r) = -i Y, (-l) k e i7rT{k+ ^e 2in(k+ ^ )z , Sr > 0, (A.l) 

k=— oo 



0i(* + 1;t) = -0i(*;t), 9 x {z + t;t) = -e~ i7rT e~ 2mz O^t). (A.2) 

We also denote 

6 2 (z-t)=9 1 {z + ^t), e A {z-T) = -ie^e i ^9 1 {z+ T --r). (A.3) 

We recall Jacobi's imaginary transformation for the theta functions: 

9l{z . T ) = _j (_ iT )-| e — 4 ^ ( - £ ; -I ) , (A.4) 

2 (z;r) = (_» r )-^e-^ 4 (- £; -I). (A.5) 

We also recall two useful summation identities (see for instance |10j): 

1 y^ e - 2mk l ei{x + y + i,r)e[{0;T) = e2mky 0i(x + ^ + fcT;nr)^(O;nr) 
n £^ 0i(x;r)0i(y + ^r) ' 0i(x + kr-nr) 0i(m/; nr) ' 

with i6Z. These two identities are equivalent through Jacobi's imaginary transformation 
(|A.4p and quasi-periodicity property (|A.lj) . 

B A determinant identity 

For two different sets of n pairwise distinct variables {x} and {y}, and 7 = Y^j=i(Vi ~ x i)' 
we consider the determinant 

det [H a ({x}, {y}) - 2Q (3 ({x}, {y})] , (B.l) 

where Qp is a n x n rank-1 matrix and 

ytt ] 1 ( gjXgj -2/j+t) gjXg - 3/j + 7 + V) ] TT gjXgj ~ Vj + ??) 

L aU > 0i (7) T 1 0i(^-%) ° 2 0i(*< -%• + */) J ^01(1// -l/i + ^) 

1 f 6>i(a?i -yj + j) 0\{xi - yj + 7 - 77) \ A Q\(xi - yj - 77) 



a 3 777 \ a 4 



(7) I " #iOj - %•) 0i (^ - j/j - »/) J jLi °i(yi - yj - v) 



^,„-„,_*V (B ' 2) 
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with a = («i, 02,03, 04) arbitrary. For t being an arbitrary complex parameter, we introduce 
the matrix Xt({x}, {y}) given by 



[X 



1 nr=i g i( x fc -yi) Gi(yj~ x k+t) Qi{x k ) 



jk Oi(t) Ui^k Oi( x k - xi) 9i(yj - x k ) Oi(x k - t) ' 
with determinant 

det [**({*}, M)]=(-l) 11 ^7-— 7Y 11^: 



We have 

det[# Q ({x},{y})- 2Q /3 ({x}, {y})] 



(*) i =1 °i( x i ~ *) f^ flifo - x fc) ' 

det n [(^ 4 F a -2^g^)({x},{y})] 



(B.3) 



(B.4) 



det n [**({*}, {y})] 



(B.5) 



We can compute the product of matrices XfH a by means of the residue theorem applied 
to the functions 



g<i> k \z) 



(z - y k + 7 + erf) A 0i(z - yi) 6i(yj - z + t) x (z) 



n 



(z-y k + efj) fj^0i(z-xi) 9\{yj - z) x (z-tY 



(B.6) 



with e € {0, +1, —1} and j, k = 1, . . . , n. These functions are doubly periodic of periods 1 and 
f, they are therefore elliptic functions and the sum of their residues inside an elementary cell 
cancels. It leads to the identities 



n n n 



£ 



n"=i 01 ( x b - J/0 01 (%' ~x b + t) X (x b ) 9 X (x 6 - y k + 7 + 677) 



fe= "t n^fc^iC^b-^) 6i(Vj- x b) 0i(x b -t) 9i(x b - y k + efj) 
h(t-y k + y + efj) ^rd^t-yi) 1 (y j ) 



i(t) 



n 



(i-M0i(t)II 



i(t-yk + tfj) iJ[ 0i(t-xi) 0i(yj -t) 

01 {yk ~y\- en) 0i {yj -y k + t + efj) X (y k - erf) 



fJl 01 (Vk -xi- efj) 0x (yj -y k + efj) X (y k -t- efj) 



nr=i0i(%- -xi) &\(vj -ty 

Gathering all terms coming from the different versions of (|B.7p . one obtains 

(X t H a )({x}, {y}) = (H a , t + 2V a , t )({x}, M)- 
Here 



(B.7) 



(Bi 



[n a , t ({x},{y})] j 



jk 



CV2 



(*) 



h {yj -yk + t + n) _ 1 

0\{yj-yk + n) 
hjyj -yk + t-fj) _ " 

0i(yj-yk-v) 



0\(yk -fj) 
h{Vk -t-fj) 

0\(yk + ff) 
h(yk -t + fj) 



<- 0\(yj) Ui^j e i{Vj ~ Vi) j j\ 9i(xi - y k + fj) yr 6i(xi - y k - fj) 

+ °i k a, i„, . _ +\ n n a. („.. _ „.\ 1 ai ll o,(. h _,„ _i_^ ° 3 11 



0i(yj- i ) IULi^i (%•-»») 



«=i 



l(^ -yk + v) 



i=i 



h(yi -yk-n) J ' 



(B.9) 
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and V a t is a rank-1 matrix: 



[v«,«(M,{i/})]. 



0' 2 



hiVk-fj) OiiVk + v) i r)l]] 

a 4— f , , ~ N I' + 0(1), 



** 2^(t)l 'hfa-t-fj) 0i(y k -t + fj) 



(B.10) 



where 0(1) stand for terms which remain finite in the limit £ — >• 0. 

Recall that t is an arbitrary parameter. We can therefore take the limit when it tends to 
to simplify the previous formula. Using the fact that V Q) i and Q^ t = ^tQa are of rank 1, 
we have 

n 

lim 0i (i) det \n a , t + 2(V a ,t ~ Qb t)] = lim 0i(i) detft a , t + 21im Vdet^JJ, (B.ll) 

6=1 

with [H® t ] jk = 9x(t)[V a ,t - Q B ,t] jk ioik = b and [H® t ] jk = [H a ,t] jk otherwise. The first 
limit in the right hand side of (|B.llj) is obviously zero, such that 

n 

lim Bx{t) det [H a ,t + 2(V a ,« - Qp, t )] = 2 Vdet [H®] 

6=1 

= det [H Q + (V Q -Qb)]~ det [ft Q - (V a - Qp)]. 
Here we have defined 



N^ = J™[^I 



jk 



, IIj/j^i(«7-w) i n 



3 jfc T-rn 



lfa -Vk + V) TT Qi(xi - Vk - V) 

° I I I 7T7 — — - «3 j 1 



U?=i°i(yj- x d I fJi 8i(yi - Vk + v) fJi Oi(yi - Vk - v) 



+ « 



i 



1(0) 



#i (%' -Vk + v) 9'i (Vj -Vk-v) 

a 27n ^\ ~ a nn ^\ 

vi [Vj -Vk + V) vi [Vj -Vk-V) 



[v a ] jk = nme 1 (t)[v a , t ] jk = 
[C/»] iik = ag«i(*)[^,t] jV 



02 — «4 



and 



W 



J k 



[Ha\jk iik^b, 

[V a -Qp) jk iffc = 6. 



(B.12) 

(B.13) 
(B.14) 

(B.15) 



Hence, 



det[fr Q ({x},{y})-2Q^({x},{y})] 



(-1)" -j-r fll (gj-Xfc) 

01(7) Wflifo-Vk) 



x { det [(H a + (V Q - Qfi))({x}, {y})] - det[(« Q - (V Q - Q/j))({x}, {y})] }. (B.16) 
If we suppose moreover that the matrix Qp is of the type 

1 



[Qb({x},M)} 



tj 



l(7) 



ft i^zi+^- A g^±i^in^±j2, ( B.i7) 



L(Si-S) 



?ite + §) J f=i *i(m + I) 
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with (3 = (/3i,/?2) arbitrary, then we can compute the product XtQp similarly as above by 
means of the residue theorem applied to the elliptic function 

J*,,_ W*-4+l) fx h(z-Vi)t>M-z + t) ^) ,r> m 



with e = ±1. It leads to 



fo+§) ft (w-§) 



fou = ss^iw t^u = a - a n t I r:7l - ( b - 19 ) 



^ 4 _>o ^ — 3» f = i^i-^i(M + ?) 
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